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ABSTRACT — The rate of muon capture in a muonic hydrogen atom is calculated 
in heavy-nucleon chiral perturbation theory up to next-to-next-to leading order. To 
this order, we present the systematic evaluation of all the corrections due to the 
QED and electroweak radiative corrections and the proton-size effect. Since the 
low-energy constants involved can be determined from other independent sources 
of information, the theory has predictive power. For the hyperfine-singlet fip capture 
rate Fq, our calculation gives Fq = 710 ± 5s~^, which is in excellent agreement 
with the experimental value obtained in a recent high-precision measurement by the 
MuCap Collaboration. 
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I. INTRODUCTION 

In a recent MuCap Collaboration experiment {l| , the rate Fq of muon capture from 
the hyperfine-singlet state of a fip atom was measured to 1 % precision. The reported 
experimental value is 

ro^P(/i-p ^ u^n) = 714.9 ± 5A{stat) ± 5.1{syst) sec"^ . (1) 

As is well known the np capture process is the primary source of information on 
the pseudoscalar form factor, Gpiq^), which appears in the nucleon matrix element 

n fl 

of the axial- vector weak current (see Eq.([6])); for recent reviews, see [3|, |J|. To be 
more specific, yup capture is sensitive to the quantity gp = G'p(g^ = — 0.88m^), where 

is the four-momentum transfer squared relevant to /ip capture (m^ is the muon 
mass). Bernard et al. used heavy baryon chiral perturbation theory (HBxPT) to 
calculate Gpiq^); their results essentially reproduce those obtained earlier by Adler 
and Dothan 6| based on PCAC, and by Wolfenstein [t] with the use of the dispersion 
relations. It is to be emphasized, however, that the systematic expansion scheme 
of HBxPT allows us to conclude that the corrections to the expression for Gp{q^) 
obtained by Bernard et al. are very small j^. The value of gp derived from HB^PT 
is gp = 8.26 ± 0.23 j^. Meanwhile, the empirical value of gp extracted from F^^'' 
with the use of the theoretical framework provided in Ref. 9| is g^'' =8.06 ± 0.55 
which is consistent with the theoretical value. 

It is to be noted that a theoretical treatment of fip capture that matches the 
1 % experimental accuracy requires a rigorous treatment of the radiative corrections 
(RCs) of order a. Czarnecki et al. ^ calculated the relevant RCs within the theo- 



retical framework developed by Sirlin and Marciano [10|, lUj • In this approach (to be 
referred to as the S-M approach), the RCs of order a are decomposed into so-called 
"outer" and "inner" corrections. The outer correction is a universal function of the 
lepton velocity and is model-independent, whereas the inner correction is affected 
by the short- distance physics and hadron structure. The inner corrections arising 
from photon and weak-boson loop diagrams are divided into high-momentum and 
low-momentum contributions. The former is evaluated in the current-quark picture, 
while the latter is estimated with the use of the phenomenological electroweak nucleon 
form factors. The expression for Fq including the RCs of order a due to Czarnecki et 



al. js] was used by the MuCap Group in deducing the above-mentioned value of 
from r'^^. We remark that, ahhough the estimates of inner corrections in the S-M 
approach are considered to be rehable to the level of accuracy quoted in the literature, 
the possibility that these estimates may contain some degree of model dependence is 
not totally excluded. This motivates us to present here a calculation of the RCs of 
order a based on model- independent effective field theory (EFT). 

In this note we evaluate the RCs for fip capture based on HB^PT, which is an 
effective low-energy field theory of QCD, see e.g., Refs. [iJQ- In HBxPT, the 
short-distance hadronic and electroweak processes are subsumed into a well-defined 
set of low-energy constants (LECs), which means that the LECs should systematically 
parametrize the inner corrections. Therefore, provided that there are sufficient sources 
of information to fix these LECs, the HBxPT approach gives model-independent re- 
sults with the possibility to estimate higher-order corrections. In two of the earlier 
publications we used the same EFT approach to evaluate RCs to order a for the 
neutrons-decay 15|, and for the inverse /3-decay reaction, i?eP e'^n, at low en- 
ergies |16,|. It is to be noted that the EFT treatments of the fip capture process, 
neutron /3-decay and the UeP — ?• e^ra reaction involve the same LECs. Therefore, if we 
can determine these LECs with the use of experimental information for one process, 
we can make model-independent predictions for observables for the other reactions. 

The remainder of this article is organized as follows. In Section II we explain the 
basic ingredients that enter into the HB^PT calculation of the fip capture rate. We 
describe in Section III the evaluation of the RCs to order a, and give in Section IV 
the numerical results for the ftp capture rate including the RCs. Finally, Section V 
is dedicated to discussion and conclusions. 



II. HBxPT CALCULATION OF THE fip CAPTURE RATE 

The theoretical framework is essentially the same as the one employed in Ref . . 
We therefore present here only a brief recapitulation of our formalism, relegating 
details to Ref. [l5]. HB^PT assumes that the characteristic four- momentum for the 
process, Q <^ ~ 1 GeV, where A^^, is the chiral scale. This theory contains two 
perturbative expansions, one in terms of the expansion parameter Q/A^ <^ 1 and 
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the other in terms of Q/itln ^ 1, where ttim is the nucleon mass. Since ttin — A^, 
the two expansions are considered simultaneously. When we include RCs in our 
considerations, a third expansion parameter a enters the theory. Our concern here is 
to carry out a HB^PT calculation up to next-to-next-to- leading order (NNLO), i.e., 
to order (Q/A^)"^ ~ a ~ 1/137. In what follows, we first describe the contributions 
that arise from the expansions in Q/A^ and Q/rriN- This part is based on the previous 



HB^PT results that can be found in Refs. 
Ref. 



17H19|. (We follow the notations used in 



17j.) We then proceed to explain our calculation of RCs. 
Muon capture being a low-energy process, the relevant weak interaction can be 
expressed as the local current-current interaction, and the transition amplitude for 
the ordinary muon capture (OMC) process in hydrogen, /x^p — > z^^ra , is given by 



In the above, Gjs = GpVud, where Gi;' is the Fermi coupling constant determined from 



the muon decay rate, and the Vud is the CKM matrix element given in Ref. [23|]. ^'^p(O) 
is the /xp-atomic wave function at r = 0. The normalization factor A/'rci, which arises 
from "matching" between the standard relativistic normalization of spinors and the 
corresponding non-relativist ic normalizations is given by A/'rci = ^T^N^Jra^Ey , where 

= = 99.149 MeV. (3) 

The non-relativistic transition amplitude, 7nr, in Eq-© is written as 

'Tnr = xlxl M Xp,Xp , (4) 

where Xp,n and x^i,l^ are the nucleon and lepton two-spinors, respectively; the explicit 
expression for the operator Ai will be given in what follows. 

The matrix element of the leptonic weak current operator, la in Eq. ([2]) is given by 
la = = Mi/7cj(l — 75)^^ , which in the present case takes the form 

^0 = ^xUl-o"- z>)Xm' ^= ^xli^-^ ■ i')^Xfi, (5) 
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where v is the unit vector in the direction of the neutrino momentum. ' 
elements of the nucleon weak current operator, j" = j"— , are given by 



'he matrix 



30| 



Jv 



W)\j':w)) 



Un{p') 
Un{p') 



m,, 



Up{p) 



(6) 



where F^^q"^), ^2(5^), and Gp{q^) are the vector, weak-magnetism, axial- 

vector and pseudoscalar form factors, respectively, and where the is the average 
nucleon mass, mjy = \{'rrtp + In the rest frame of the initial proton, the non- 

relativistic nucleon currents in HBxPT are given by[31l] 

2m7v 



3: = Mnn,{p') 
1 



E'+rriN 



X 



+ 



E'+niN 
2 



%^^F;{q') 
E'+rriM ^ ^ . 

1 



E' + rriN 
j: = ArnMp')lGA{q') 



2mN 

[S^,S-q]{F^q^)+F^,{q'))\pM 
2{S-q) 



2S, 



G 



E' + niN 

2. 2{S-q)qa 



m^{E'+mM) 



PviO) 



with the heavy nucleon spinors, ny{r') and Pt,(0) defined as jl2| 

1 



nv{p') 



lii+i')Mp') , pM = . 



(7) 



(8) 



The kinematics in the rest-frame of the proton is as follows. The four-momenta 
of the initial proton and the outgoing neutron are p = {mN,0) and p' = {E',p'), 
respectively, where E'= ^/m^ + p'^ and p' = —pu- The four-momentum transfer in 
the OMC process is g = p' —p = {qo,c\), with go = E' — miy = + 0{mjj^) and 
q = —pi, . Expanding the proton and neutron spinors in Eq.(I7]) up to 0{m~^) leads 
to 

j(g) = -xi W X ^)/2 (g) + ^/s (g) + c?/r(g) + v{a ■ u)f^{q)] x, , (9) 

where the non-relativistic polar- vector form factors are related to the standard Lorentz 
covariant form factors in the proton rest frame via 



fiiq) = F^iq'' 



q' 



im 



N 



+ 



F^iq'' 



4m% 



(10) 



while the non-relativistic axial-vector form factors are related to the covariant axial 
form factors via 



N 



f2{q) 



|q| 



1 + 



5m 



N 



Gp{q'). (11) 



nil) 



GA{q') + 



Gp{q' 



(12) 



2mN V 2m^m]\f 
The non-relativistic form factors appearing in Eqs. (ITOl) . ( ITTl) and ( |T2l) have been 
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2l| . up to next-to-next-to leading order (NNLO) or 



calculated in Refs. 
0{{Q/A^)^), in HBxPT. In the proton rest-frame, they are given by 

,2 



q 



4m% 



(4vr/.) 



qHl9\ + 2BrA + qHl9\ + -]\n 



dz 







Ml{^g\ + l) - q^z{l-z){bgl + l) 



In 



l-zil-z) 



Ml 



|q| 



2m N 

|q| 

2m AT 



/2(?) 
/3(?) 

/r(g)=^?A(i 

/2(9) 
/3(g) 



1 + + fl-A 



(4vr/, 



1 -zfl-zl 



Ml 



N 



+ 



|q| 



|q' 



(47r/.)2 

,2 



B 



3 ; 



9 



5m 



2M2 



(47rA 



r52 + 



|q| 



(4vr^) 



ri^3, 



2m 



1 



TV 



(13) 



In terms of the quantities derived above, the operator [see Eq.([2])] is written as 



M = {l-ayO)[fl{q)-iar{BN^y)r2{q)-{^ry)n{q) 

-{^i-^N)n{q) - {^i ■ ^){ffN-i^)f2{q) + {^N-v)n{q). 



(14) 



where Bi and (Tn are the spin matrices acting on the lepton and nucleon spinors, 
respectively. We may choose the direction of the emitted neutrino as our 2;-axis, i.e.. 



6 



1/ = z. In the helicity basis, the amphtude 7nr appearing in Eqs.(l2]) and ^ is given 
as 

six 

where h = L {h = R) corresponds to the positive (negative) hehcity state of the final- 
state neutron, and S = (S = 1) represents the hyperfine-singlet (triplet) state of the 



muonic hydrogen atom [32l]. The constraint, = S"" — |, reduces the eight possible 
helicity amplitudes in Eq. ffTSj) to the following three: 

M(L;0,0) = V2{r,+2f^ + r3+3f^ + f^ + f^), 
M(L; 1,0) = V2{r,-2f^ + r,-f^ + f^ + f^), 

M{R; 1, -1) = 2 (/[ + r3 - n - n - n) ■ (le) 

Finally, since the binding energy of the muonic hydrogen atom can safely be ignored, 
the total OMC rate in a hyperfine state S is given as 



2{m,+m,) 2Sri j {2.Y2E' [2.f2E^ ^'"^'^'^"^^ ' ^' ' l"^^-' 

S z jh 

2 / \ 

Y,\M{h;S,S^)\ (17) 



2S + 1 GAnm^mN \E^+^m% + E^ 



Sz,h 



where Pj is the initial total four-momentum. If we ignore radiative corrections and 
identify "^^pi^O) with the lowest-order Is-state Coulomb wave function, $is(0) = 



a^U'^p/TT)"'" with iXjip = m^,mp/ (m,,+mp), then the last line in Eq. (fT7|) agrees with 



^'l^lp/^V' wit 

Eq.(26) in Ref. 17|]. In the present work, however, we do include radiative correc- 
tions, and it turns out that, at 0{a) under consideration, there appear two types of 
significant radiative corrections to $1^(0), and these corrections will be discussed in 
the following section. 

We now evaluate the form factors in Eq. flT^ . which determine the helicity am- 
plitudes in Eq. fll6l) . Table [T] shows the numerical values of the nucleon weak 
form factors and helicity amplitudes calculated for the four-momentum transfer, 

= ql = — 0.88m^, relevant to OMC. These numerical values were obtained with 



the use of the following input parameters: gA = 1.266, Ky = 3.706, /tt = 92.42 MeV, 
= 139.57 MeV, and m/^=_938.919 MeV. The LECs appearing in Eq.^ are 



determined following Refs. 
(G-T) discrepancy relation. 
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20|. First, i?2 is fixed from the Goldberger-Treiman 



2M^ 



Bo 



gA rriN 

g-nNNf-K 



1. 



:i8) 



For g^NN = 13.40 and gA = 1.266 (see, e.g., PDG2002 |2^), this relation leads to 
i?2 = —1.761. The values of gA and g-j^NN have been slightly changing over the years; 
if we use the latest values g^NN = 13.05 and gA = 1.270 (taken from PDG2012 23l|). 
we obtain B2 = —0.498. To what extent the existing uncertainties in qa and gTrNN 
affect the calculated fxp capture rate will be discussed in the last section. The LEG, 
i?3, is fixed from the nucleon axial radius. 



The value of the iso- vector axial radius {r\)^^'^ has large uncertainty, see e.g., Ref. 12|. 
From the empirical axial form factor G'yi(t) = (7^1/(1 — t/m^)^, we find = 
"U/rriA = 0.62 fm (0.57 fm) for mA = 1100 MeV (1200 MeV). We adopt the value 



(r2 )i/2 = 0.65 fm cited in Ref. |2oj to find B3 = 3.08. The last of the LEGs in Eq.(fT3|). 



b[^q , is related to the nucleon iso- vector form factor 



20| 



(471^)2 6(47rA) 



5^7i + l m 
3(4vrA)2 



From the measured value of (r^) we obtain B^''' 
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0.27 for Av = 1 GeV. These values 



of the LEGs were used in obtaining the numerical results given in Table HI 



TABLE I: The OMC form factors Eq.([T3]) and helicity amplitudes Eq.([l6]) at ql 
obtained for qa = 1.266 and Qt^nn = 13.40. 



-0.88m2, 



/r(g*) 










nil*) 


M(L;0,0) 


M(L; 1,0) 


M(L;1,-1) 


0.966 


0.244 


0.053 


1.245 


-0.419 


0.044 


3.447 


-0.770 


0.148 
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III. RADIATIVE CORRECTIONS 



In this section we consider radiative corrections to OMC, which consist of the 
usual QED loop corrections and loop corrections involving a weak-interaction vertex. 
Relegating the discussion of the latter to the end of the section, we first discuss the 
QED loop corrections. 

The initial state in fi~p capture is a charge-neutral yU^p-atom, and the final state 
involves only electrically neutral particles. Therefore, to the order in HB^PT under 
consideration, the bremsstrahlung process does not contribute to the "standard" ra- 
diative corrections (we ignore the higher order, i.e., 0{a/mi<^) corrections which are 
negligible.) There are, however, two QED loop corrections to the initial state wave 
function which must be considered: the vacuum polarization correction, Sip^FW^ 
the correction due to the finite proton size, (5'?/'fjf (0). Inclusion of these corrections 
changes the lowest-order muonic atomic wave function, $1^(0), into \I'is(0): 



(19) 



Eiras and Soto 



24j calculated SipY^{0) to order 0{a), while Friar [25|] discussed 0{a) 
contributions to (J'j/'ff (0). 

n 

The analytic expression for 6iIji,^{0) derived by Eiras and Soto [2^] reads 



W(0) = - 



+ 



a 

IT 
11 



^-^e + ^^^-^e^ + ^(e^+e^-2)F,(0 



i 



2tt 



6 



6(^2-1) 



TT 



ISvr 



+ ^3 + 4^-9^ + 1^ 



-(i3e'-iie'-ii)Fi(o 



2, in,, a 



multi— 7 



pole 



(20) 



where ^ = me/{afi^p) ~C(1); the expressions for the functions -Fi(0 (^ = 1)2,3) in 



Eq. fl20|) can be found in Ref. 



24| . As explained in Ref. 2J], the first curly bracket 



in Eq. ( l20l) corresponds to zero photon exchange contributions, the second bracket 
corresponds to Coulomb pole subtraction terms, and the third bracket represents the 
multi-photon exchange contributions. Thus 6tpYFi^) consists of three parts: 



Si^l^iO) = [Si^l^'m^ + [<5^r/(0)]pole + [5^r(0)]n.ulti 



(21) 



We denote by the fip capture rate for the hyperfine-state S {S = or 1) , obtained 
by using $i,(0) for ^^,^(0) in Eq.(^. The use of $i,(0) [1 + 5^^/(0)] for ^^,^(0) in 
Eq. f|T7|) changes r?-* into 



(22) 



Table Hn shows the numerical results for ((5r5)^^/r|'-* = 26tljYF ■ '^^^ ^^^^ three columns 



show the individual contributions of the three terms in Eg. (121]) . while the fourth 
column gives 25-?/'^^, which is the sum of these three contributions. For comparison, 
in the fifth and sixth columns, we quote the values of 26tljY.F (i^ notation) obtained 
m Refs. Igl,^- 

Our result for 2 6iIjY^{0) agrees with the value given by Czarnecki 



TABLE II: Corrections from vacuum polarization (VP) effects, ((5r,s)^^/r 



The last two columns give the values of 25^7/(0) in Refs. [9|, 



2#vP(o 



261 ] for comparison. 



Zero photon 


Coulomb pole 


Multi-photon 


Total VP 


Czarnecki, 




exchange 


subtraction 


exchange 


contribution 


Marciano, 


Goldman 


2^7/(0)]^ 


2W(0)]pole 


2[<5<P(0)]„,,iti_^ 


2 5<P(0) 


& Sirlin [9] 


[26] 


1.045^ 

7r 


0.358^ 

n 


0.250^^ 

TT 


1.654^^ 

7r 


TT 


2.95^ 

TT 



et al. jol within ~5%. Since the size of the 26ipYF correction itself is about 0.4%, 
we can say this part of QED corrections is controlled with sufficient accuracy for our 
purpose. 



The proton finite-size correction up to 0{a) is given as 25] 



5<^(0) = -a/i^p(r)p. (23) 

where (r)p is the first moment of the proton charge distribution, Pp(r). Unfortunately, 
{r)p cannot be measured directly, whereas the second moment, (r^)p, can be extracted 
from experimental data. In order to evaluate (r)p, we assume a certain functional 
form of the proton charge distribution, Pp(r), involving a single parameter, and after 
determining this parameter from the measured value of (r^)p, we deduce (r)p from the 
assumed Pp(r). Table UTTl gives (r)p and (r^)p calculated for three different functional 
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forms of Pp(r). The results for the exponential form, Pp(r) = 1/ (87rrQ)e~'^'" are given 
in the fourth column; the exponential form corresponds to a dipole-type proton form 
factor (in momentum space), which reproduces very well the elastic electron-proton 
scattering data. We also present the results for two other commonly used forms 
for /5p(r), the uniform distribution (second column), and the Gaussian form (third 
column); these results have been extracted from Ref. {25 1. The last row in table UTTl 



TABLE III: First and second moments of the proton charge distribution calculated for 
various forms of Pp{r) characterized by a single parameter. 



Pp(r) 


Uniform 


Gaussian 

( 1 ]\-{r/rof 


Exponential 


{r)p 


3R/4: 


3rl/2 


3ro 
12 




VTb/i = 0.968 


2y^2/37r = 0.931 


V3/2 = 0.866 



shows the ratio (vp) /^/{r^p for each assumed form of Pp{r). By taking the average of 
the results for these three cases, we deduce {r)p = (0.916 ± 0.051) /p) ; the ' error 
estimate" here has been obtained by interpreting the scatter of the results in table IIIII 
as a measure of uncertainty. Then, with the use of the experimental value of the 
proton r.m.s. radius, VT^p = 0.862 fm we find (r)p = 0.790 ± 0.044 fm. 

Using this value in Eq.(l23]) leads to (0) ~ -0.00275(1 ±0.056). Correspondingly, 
the finite-proton-size correction to the capture rate in Eq. (fT7|) is found to be 
2 6tp^f{0) = —0.0055(1 ± 0.06). This result is essentially the same as that given in 
Eq.(8) of Ref. [9|]. Thus, the finite-proton-size correction is of the same order as the 
vacuum polarization correction shown in Table [Tll 

In addition to the two QED corrections discussed above, we need to consider the 
"standard" radiative corrections involving a weak-interaction vertex. It is to be noted 
that part of these corrections are already included in Gp, if one uses (as we do here) 
the value of Gp determined from the measured muon lifetime. In the following, what 
we simply call the "electroweak loop corrections" refer to those electroweak loop 



11 



corrections that have not been accounted for by the use of the Gp derived from the 
muon hfetime. We remark that, to the order in HBxPT under consideration, the 
electroweak loop corrections are identical for np capture and neutron beta-decay. We 



can therefore utilize the results obtained for neutron beta-decay in, e.g., Refs. |15l.l29|. 
Since the muon velocity, /3, in the initial /ip-atomic state is essentially zero, we can 
take the limit of /3 — )■ in t he p revious evaluations of the radiative corrections to 



the neutron beta-decay rate 



15 



291], (In applying the results obtained for neutron 



/3-decay to the /ip capture case, we must drop the bremsstrahlung contributions, since 
both the initial and final states in /ip capture contain only charge- neutral particles.) 
Then the electroweak radiative loop correction to the [ip capture rate is obtained as 

r^^ — )• r^^(i + rCew) 1 



with 



RG 



EW 



a 
2^ 



evim^) +3 In 



TUN 


271 







(24) 
(25) 



In this expression the electroweak LEG, ey(m7v), subsumes short-distance physics not 
probed in the low-energy muon capture reaction. The value of this LEG at the scale, 
A = rriTv, has been determined in Refs. 15|, ll6|] by comparing with the expressions 
for the short-distance radiative corrections derived by Sirlin and Marciano 10|, 111] 
for the electroweak processes. The result is eyirriN) = 19.5. In the next section we 
discuss the numerical consequences of our evaluation of the 0{a) radiative and finite 
proton-size corrections discussed in this section. 



IV. NUMERICAL RESULTS FOR THE CAPTURE RATES, Tq AND Ti 

As explained earlier, Fq^'' (vf'^) denotes the hyperfine-singlet (hyperfine-triplet) 
capture rate calculated without including radiative corrections; viz., Fq^^ and f|^°^ are 
obtained by identifying \E'^p(0) in Eq.l fTTI) with $is(0). Using the inputs listed in 
Table [H we obtain 

f[,°^ = 693 s"^ and vf^ = 12.0 s'S (26) 

corresponding to the use of = 1.266 and gT,NN = 13.40. The inclusion of the radiative 
corrections discussed in Section III modifies F^^ {S = 0, 1) to F5 as 

Ts = rfil + i?CQED + RCew) ; S = 0, 1. (27) 

12 



Here RCqed represents the corrections arising from the change in the atomic /ip 
wave function due to the vacuum-polarization and finite-proton-size effects, while, as 
explained earlier, RCew is the electroweak radiative correction: 



RC^EB = 2S^P^F{0) + 25^ff (0) (28) 



a 



RCew = —\4imN) + S\n 



'ruN 


271 







(29) 



We remark that, since the last two terms in Eg. (1291) almost cancel each other, RCew 
has a pronounced dependence on the LEG, ey{mN), which characterizes the short- 
distance processes. 

The numerical consequences of including the radiative corrections are displayed 
in Table HV] where Ts (5 = 0, 1) are shown along with and the changes due to 
the individual contributions of RCqed and RCew- Again, these results have been 
obtained with the use of qa = 1.266 and gnNN = 13.40. Table [IV] demonstrates that 
the largest radiative correction to the OMC rate comes from RCew, in conformity 
with the results reported in Ref. jsj. In particular, for the hyperfine-singlet OMC 
rate, which is of our main concern, RCew changes Fq^'' by ~2 %. 

TABLE IV: The hyperfine-singlet and -triplet OMC rates, Tq and Ti (in units of s~^), 
calculated with and without radiative corrections (the proton-finite-size effect is included 
as part of RCqed) corresponding to ^TrTVAf = 13.40 and 3^ = 1-266. 



p(0) 
^ 


ri°\i+i?CQED) 


T^^\i + RCew) 


To = ^o*\'^+RCqed+RCew) 


692.9 


691.7 


708.4 


707.2 


r(o) 
J- 1 


Tf\l + RCQED) 


T^Xi + RCew) 


Ti = Tf\l+RCQED+RCEw) 


12.0 


11.9 


12.2 


12.2 
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V. DISCUSSION AND CONCLUSIONS 



In the previous section we have presented our numerical results obtained with the 
use of representative values for the relevant input parameters. We now discuss to what 
extent the uncertainties in these input parameters affect the calculated values of the 
/ip capture rates, Ts {S = 0, 1). We shall chiefly concentrate on the hyperfine-singlet 
rate Tq, a quantity of primary concern for most fip capture experiments. 

As mentioned, the LEG, B2, is determined from the G-T discrepancy [see Eq.( [T8|) ]. 
and the fact that the current precision of the values of qa and Qt^nn is somewhat 
limited leads to rather significant uncertainty in B2. The results in Table HVl were 
obtained for B2 = —1.761, which corresponds to Qt^^n = 13.40 and (7^ = 1-266 taken 
from PDG2002 {22]. If we adopt Qt^nn = 13.05 and gA = 1.270 (values given in 
PDG2012 [23]), then we obtain B2 = -0.498 and, correspondingly, r[,°^ = 699 s"^ 
and r[*''* = 11.7 s~^. Thus, the uncertainty in B2 changes by 7 (about 1% 
increase), and vf'^ by 0.2 s"^ (about 2 % decrease). If we take into account (in the 
last column in Table WV} the mentioned variation in rQ^"*, the corresponding change 
in Tq ranges from 707.2 s^^ to 713.7 s^^; thus 

To = 710.4x(l± 0.005) (30) 

where the relative error was deduced from the 1 % difference between the above- 
quoted two values of Fg^-* . 

We next consider the uncertainty in the proton axial radius, ((^'i))^''^, discussed 
in Section II. The results shown in Table HVl were obtained for ((r^))^^^ = 0.65 fm. 
If we instead use ((r^))^^^ = 0.57 fm, corresponding to = 1200 MeV, we find 
ro^"* = 695. 7s~^, an increase of 2.9 (or ~ 0.5%). Again, if we consider (in the last 
column in Table [IV]) the scatter in the value of T^q \ then the corresponding change 
in Fq ranges from 707.2 to 710.1 s~^, i.e., 

Fo = 708.6 x(l ±0.0025)5-^ (31) 

where the relative error was deduced from the 0.5 % variation in Fq^''. Taking the 
average of the values in Eqs.( l30|) and ( 13T1) . we arrive at 

Fo = 710x(l± 0.006) (32) 
14 



where the error has been deduced from the quadratic sum of the errors in Eqs.f l30p 
and (|3T|). 

In connection with Table [IV] we have pointed out that, of all the corrections of 
0{a), the electroweak loop correction, RCew^ is largest; it increases V'^^ by as much 
as ~2 %. So, if RCew is not evaluated with sufficient accuracy, the theoretical error 
in To can be larger than indicated by Eq. (l32!) . As already mentioned, RCew is a 
sensitive function of the LEG, ey(Mjv), as the last two terms in Eq. fl29l) nearly cancel 
each other. In the present work, following Ref.jl^, we have determined eyiM^) by 
comparing our HBxPT results with those obtained in the S-M method 10|, • Since 



N) 



this method is generally considered to be highly reliable, we believe that ey{M^ 
is known with sufficient accuracy to make the uncertainty in Fq related to RCew 
much smaller than 0.6 %, the error arising from the other sources [see E g. (321) ]. We 



remark that the same LEG, ey{mj^), also appears in neutron beta decay [15|| and the 



inverse beta decay process, Ue+P — ^ e~^+n [16|. It is therefore, in principle, possible 
to use either the neutron /3-decay or the fip capture to control ey(mjv) and make 
predictions for the other processes involving the same LEG. This would allow us to 
deduce e(?(mAr) without using the result of the S-M method. 

In conclusion, the present HB^PT calculation of the hyperfine-singlet capture 
rate Fq, including radiative and proton finite-size corrections of 0{a), gives 

Fo = 710±5s-^ (33) 

This is in excellent agreement with the experimental value quoted in Eq.([T]). The 
0.6 % theoretical error in Eq. (l33|) is dominated by the uncertainties in the input 
values of qa and gnNN that enter into the G-T discrepancy. 
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